Diffusivity and configurational entropy maxima in short range attractive colloids 
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We study tagged particle diffusion at large packing fractions, for a model of particles interacting 
with a generalized Lennard- Jones 2n-n potential, with large n. The resulting short-range potential 
mimics interactions in colloidal systems. In agreement with previous calculations for short-range 
potential, we observe a diffusivity maximum as a function of temperature. By studying the temper- 
ature dependence of the configurational entropy — which we evaluate with two different methods 
— we show that a configurational entropy maximum is observed at a temperature close to that of 
the diffusivity maximum. Our findings suggest a relation between dynamics and number of distinct 
states for short-range potentials. 

PACS numbers: 82.70.Dd, 64.70.Pf, 61.20.-p 
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In recent years, several studies have been focused on 
short-range attractive systems, in which the range of the 
attractive part of the potential is much shorter then the 
typical repulsive core. In nature this behavior is en- 
countered in colloidal systems, in which the range of the 
attraction can be finely tuned by changing the solvent 
properties or the chemistry of the dispersed particles. 
The short length-scale of the attraction generates many 
unusualphenomena, (for recent reviews see for example 
From a thermodynamical point of view, it is 
now well established that the liquid-liquid coexistence be- 
comes metastable with respect to the fluid-solid one (^(^ , 
a phenomenon which is not observed in atomic or molec- 
ular systems where the interaction is always long-ranged. 
Recently, theoretical, numerical and experimental stud- 
ies have analyzed the dynamical properties of these sys- 
tems. One of the most astonishing discoveries is that, at 
high density, the metastable liquid is characterized by a 
non-monotonic temperature (T) dependence of the diffu- 
sivity: dynamics slows down not only upon cooling (as 
is commonly observed in molecular systems), but also 
upon heating. The slowing down upon heating can be so 
intense that a novel mechanism of arrest takes place at 
high T , creating in the same system two different glass 
phases: one at low T called the attractive glass and one 
at high T, the repulsive glass. The former is generated 
by the short-range attractive part of the potential, while 
the latter is due to the cage effect, i.e., caging of the par- 
ticle in a shell of neighbors, analogous to the hard-sphere 
lass. This scenario has been confirmed experimentally 
I IM IS and by numerical simulations [iCt . 

The existence of a fluid phase between two glass re- 
gions suggests that the diffusion coefficient D exhibits 
a maximum when the strength of the attraction (or the 
temperature) is varied along a path of constant density. 
Such a feature has been observed in simulations of short- 
range square well models. The non-monotonic behavior 
of D has been explained in term of Mode Coupling The- 
ory as resulting from the competition of two different 



localization lengths, one associated with the "hard-core" 
and one with the short-range attractive bond fl5|. In 
this Letter we calculate the T dependence of the config- 
urational entropy Sconf — a measure of the number of 
distinct states of the system — with the aim of provid- 
ing insights into the physical origin of the D maximum 
■ We use two different routes to determine Sconf '■ the 
first one, based on a Potential Energy Landscape (PEL) 
investigation, requires an estimate of the vibrational free 
energy of the system close to the explored local minima 
of the PEL (the so-called inherent structures (IS) 14] ); 
the second one uses a perturbed Hamiltonian 15, 16] and 
requires the mean-square distances from a given equilib- 
rium configuration in the limit of vanishing perturbation. 
In both cases, the configurational entropy is calculated 
as a difference between total entropy (estimated via ther- 
modynamic integration from the ideal gas state) and the 
vibrational entropy. In the PEL formalism, it is crucial 
that the interaction potential of the system is continuous, 
so that IS can be properly located via a steepest descent 
minimization of the potential energy and so that vibra- 
tional frequencies can be properly calculated. For this 
reason we focus on a continuous model that possesses a 
steep repulsion and a short-ranged attraction, and that 
has been proved to reproduce features of short-ranged at- 
tractive colloidal system discussed above [T^. The main 
results of the present work are: i) the presence of a diffu- 
sivity maximum on varying the temperature for the con- 
sidered attractive-colloid model that confirms that origin 
of this phenomenon can be ascribed to the potential being 
short range rather then on the particular shape chosen 
|lOl ; ii) the existence of a maximum in the T depen- 
dence of Sconf — independent of the chosen method to 
estimate it; Hi) the fact that the temperature at which 
D and Sconf have a maximum nearly coincide. 



The chosen model of attractive colloid is based on a 
generalization of the Lennard- Jones pair potential (LJ 
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2n-n) proposed by Vliegenthart et al. 



takes into account the number of distinct local states 
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(1) 



^(r, p) = Sconf{T, p) + Si,ib{T, p) 



(2) 



Differently from the original proposition we choose an 
extremely large value of the exponents, i.e. n — 100. For 
this value of n, the potential (reported in the inset of 
Fig. ^ provides a continuous version of the 3% square- 
well potential, which has been extensively studied in re- 
cent years 

We perform standard isothermal molecular dynamics 
using the Nose-Hoover thermostat JSJ. The simulated 
system is composed of 256 particle confined in a cubic box 
with periodic boundary conditions. In order to prevent 
crystallization we study a 50 : 50 binary mixture with 
the following parameters: aAA/o's_B=1.2, <jab — {o'aa + 
o'bb)I'2-, ^AA=eBB=£AB- We use Lennard Jones units 
{aAA for length, caa for energy, T^{ma\j^/eAAV^'^ for 
time). We chose the Boltzmann constant ks — I, conse- 
quently the temperature is measured in unit of caa- Cut 
and shift in the pair potential are used (rcut= 1.4). The 
time step used is ^0 = 7x10"^. The packing fraction in- 
vestigated is 0=0.59, corresponding to density p ~ 1.43. 

Fig. ^ shows the mean-square displacement, 
r2(t)=Ar-i([r(t) - r(0)]2), where r is the 37V vec- 
tor of the coordinates, for different temperatures ranging 
from T=0.26 to T=2.0. Result are in agreement with 
previous findings based on the square- well model [lH : 
at high T there is a well defined plateau at about 
X 10^^, while at low T the plateau disappears 
and one observes a transient sub-diffusive regime before 
the diffusive one. It is worth noting that the value of 

~ 10""^ at the beginning of the transient regime at 
low T corresponds to the length of the attractive well 
of the pair potential. Fig. shows the D (evaluated 
from the long time limit behavior of and rescaled by 
the quantity Da=^J a'^T /m) as a function of T. One 
observes a maximum located at about Tmax = 0.4. 
On lowering T, for T < T,nax, D decreases quickly, as 
the system approaches the attractive glass line (almost 
horizontal in the cj) — T plane). On increasing T, for 
T > T,nax, D decreases smoothly, consistent with the 
observation that the repulsive glass line is almost vertical 
in the <f> — T phase diagram. 

In analogy with recent studies for atomic and molecu- 
lar liquids, one can ask if the T dependence of D close to 
arrested states is correlated to the T dependence of the 
Sconf, as, for example, in the well known Adam-Gibbs 
proposition In the case of short-range models, a 

test of the correlation between D and Sconf can be ex- 
ploited in a more direct way, capitalizing on the presence 
of the _D-extremum at T = T^ax- We try to answer this 
question, calculating the configurational entropy for our 
model system. To estimate Sconf, we write the total en- 
tropy as a sum of two contributions: a local vibration 
entropy Syit and a configurational entropy Sconf that 



This expression, which can be formally derived within the 
PEL framework |lj| and within mean field for models 
of disordered p-spin systems J9j, is based on the idea 
that dynamics is described by two well separated time- 
scales: a fast dynamics describing local rearrangements of 
particles within a state (vibrations around local minima 
of the PEL) and a slow dynamics which accounts for the 
slow exploration of different states. 

Eq. 12 shows that Sconf can be calculated from the 
knowledge of total entropy S and vibrational entropy 
Svib- The total entropy S can be calculated using thermo- 
dynamic integration along paths in the T-p plane. With- 
out entering in details (see Ref. [13), we can write, mea- 
suring entropy in units of fc^: 



(3) 



where U is the potential energy, Tq is a reference temper- 
ature (ro=0.4 in our case), and S{To,p) is given by the 
following equation, which describes the path at constant 
Tq connecting the ideal gas to the reference state (Tq, p), 



SiTo,p)^S'^{To,p) + 
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where S'*'' is the ideal gas contribution (which includes 
the entropy of mixing, since we are dealing with a binary 
mixture), and Pex is the excess pressure. Performing nu- 
merical simulations at 15 different p values from 0.036 up 
to 1.434, we estimated Pex{To, p) and U {Tq, p), which, to- 
gether with U (T, p) allow us to calculate with sufficiently 
high precision the T dependence of S. The calculation 
of the vibrational entropy requires some approximations. 
For this reason we follow two independent different cal- 
culation routes. In the PEL approach, the vibrational 
entropy can be calculated from the local curvature of the 
PEL around the explored inherent structures In 
harmonic approximation, the vibrational entropy can be 
written as (omitting the density dependence): 



SilliT)^ [l-ln(/3fi^.)] 



(5) 



where h is the Plank's constant, f3 — l/fc^T and tOi are 
the eigenfrequencies of the Hessian matrix evaluated at 
the inherent structure. The T dependence, besides the 
factor /3, is contained in the different properties of the 
local curvatures around the inherent structures explored 
at different temperatures (i.e. in the T dependence of the 
uji values) lU]. 

Fig. 1213 shows the configurational entropy per parti- 



cle S^^l^^/N, calculated as the difference between 



total 
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entropy in Eq. 12 and vibrational entropy S'^]l given by 
Eq. The T dependence shows a maximum at tem- 
perature Tmlx — 0.6, shghtly higher than that of dif- 
fusivity {T„iax—0-'i)- The presence of the maximum is 
a remarkable result, indicating a close relationship be- 
tween temperature dependence of diffusivity and config- 
urational entropy, even if a quantitative coincidence of 
the peaks seems not to be exactly achieved. Of course, 
the use of the harmonic approximation deserves a few re- 
marks: while in the simple liquid the harmonic approxi- 
mation works well at sufhciently low T, in our case, due 
to the steepness of the pair potential, the harmonic ap- 
proximation is expected to break down at a much lower 
T. Moreover, an estimate of the anharmonic contribu- 
tions, following the techniques developed for atomic and 
molecular systems is not feasible in the present case, due 
to the strong T dependence of the anharmonic energy. 
In order to corroborate the above finding, we calculate 
Svib using an alternative method, based on a Perturbed 
Hamiltonian approach [isl Il6l l22l| . One considers a per- 
turbed Hamiltonian: 

pH' = (3H + a{r - rof , (6) 

where H is the original Hamiltonian, a is the strength of 
the perturbation and Tq is a given equilibrium reference 
configuration. The free energies F{a) of two systems 
with different a values (ctoo and aa) are related by: 

pFiaoo) = f3F{ao) + j ^ da' {{r - Vof)^, , (7) 

J CIO 

where (...)«' is the canonical average for a specified a'. In 
the large aoo oo limit, F{aoo) has the exact expression 

3iV 

l3F{aoo) = In X + l3Eo + — ln(aoo/7r) , (8) 

where Eq is the potential energy of the reference config- 
uration ro and A is the thermal deBroglie wavelength 
A = (27r/3S^/m)^/^. If a small ao value can be cho- 
sen in such a way that the corresponding system is 
equivalent to the original system, but constrained to ex- 
plore only the phase space of one state, the estimate of 
da' {{r — ro)^)a' is sufficient to evaluate the required 
vibrational free energy. 

Writing the free energy as a sum of a potential energy 
term and an entropic term, (]F{ao) = (3Eq + 3N/2 — 
Syib, the following expression for the vibrational entropy 
is derived: 

S^(T)- /7^"'«r-rof).'-^ln(^^) , (9) 

where e is the Neper number. We use isothermal molec- 
ular dynamics with Hamiltonian H' to calculate ((r — 
ro)^)Q at different a and T. We perform averages over 



20 different reference configurations rg, chosen from equi- 
librated configurations with unperturbed Hamiltonian H 
at temperature T. Fig.Elshows ((r — ro)^)Q as a function 
of a for different T. The dashed line is the (T indepen- 
dent) high a limit, 3/(2a). As discussed above, one has 
to choose an ao value in Eq. |3 in such a way that the 
system remains trapped in a given local state. While 
at high temperature the data are smooth and the values 
of ((r — ro)^)a remain well below the cage value (about 
~ 2 X 10~^), at low T the behavior is quite different: 
starting from high a, we observe first an approach to a 
small value r'^ 2 x 10^"^ (corresponding to more or less 
the same behavior in the mean square displacement, even 
if less pronounced - see Fig. ^ and then a departure from 
it at smaller values of a. We interpret the former as the 
vibrational motion inside the state. The full line in Fig.|3| 
is a guide to the eyes that extrapolates the first behavior 
at lower a for the T=0.26 case. We have then chosen as 
ao in Eq. El a value for which the system has not yet left 
the Hue: ao = 5x 10^ (indicated by an arrow in Fig. ISJ. 
Although this is a feature only of the low T data, we have 

(2) 

chosen the same ao for the estimation of S^^f^ for all T, in 
order to obtain a coherent definition of it. The aoo has 
been fixed at 2 x 10^, where ((r — ro)^)Q has reached the 
asymptotic behavior (dashed line in Fig. 13)1. Fig.Et shows 
the T dependence of the configurational entropy per par- 
tide S^J^j:/N, calculated as the difference between total 
entropy and S^^l . Again one observes a peak, located 

(2) 

at about Tmax — 0.5, close to that of diffusivity. We note 
that using the same value of ao for all the temperatures 

(2) 

introduces an underestimation of S^^j^ , more pronounced 
for the high T data. This fact could have the effect of 
moving the peak to a lower T value, approaching the 
peak value of the dilfusivity. Although the absolute val- 
ues of S'^o^y and s'^^j: are different, the T dependence is 
quite similar, suggesting that the errors are mostly T in- 
dependent. More interestingly, in both approximations, 
a maximum in Sconf , not far from the T at which D has a 
maximum, is clearly detected. Since both methods give a 
T dependence of the configurational entropy with a peak 
located close to that of diffusivity, our work strongly sup- 
port the possibility that in short-range colloidal systems 
the diffusivity maximum is related to a maximum in the 
number of states visited by the system. 
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FIG. 1: Mean square displacements for different temperatures 
(from T'=0.26 to T=2.0). In the inset the pair potential LJ 
2n-n, with n=100. 
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FIG. 2: a) Diffusivity D as a function of temperature. The 
value of D is plotted rescaled by Dq = \J a'^T jm. b) Config- 
urational entropy 5^^^^ per particle using Eq. |K1 to estimate 
the vibrational entropy {PEL approach), c) S^^j per parti- 
cle using Eqjpior vibrational entropy {Perturbed Hamiltonian 
approach) |2^. 
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